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ABSTRACT 

We examined X-ray spectral and timing properties of ultraluminous X-ray sources (ULXs) in nearby 
galaxies in XMM-Newton archival data. There appear to be three distinct classes of spectra. One 
class shows emission from hot, diffuse plasma. This thermal emission is similar to that seen from 
recent supernovae; the temperatures are in the range 0.6-0.8 keV and the luminosities are the lowest 
in our sample, near 10'^^ erg/s. Three sources have spectra which are strongly curved at high energies 
and have the highest temperatures in our sample, 1.0-1.4 keV. These spectra are well fitted with a 
power-law plus multicolor disk blackbody model with the power-law dominant at low energies or a 
Comptonization model. The remainder of the sources are best fitted with a power-law plus multicolor 
disk blackbody model, as is commonly used to describe the spectra of accreting black holes. These 
sources have the lowest thermal component temperatures, 0.1-0.4 keV, and extend to the highest 
luminosities, above 10^" erg/s. The temperature of the thermal component is in three distinct ranges 
for the three source classes. This diversity of spectral shapes and the fact that the sources lie in 
three distinct temperature ranges suggests that the ULXs are a diverse population. Two ULXs which 
show state transitions stay within a single class over the course of the transition. However, we cannot 
conclude with certainty that the classes represent distinct types of objects rather than spectral states 
of a single population of objects. More monitoring observations of ULXs with XMM-Newton are 
required. We also searched for timing noise from the sources and report detection of noise above the 
Poisson level from five sources. In three of the sources, the power density spectrum increases with 
decreasing frequency as a power-law down to the lowest frequencies observed, below lO^** Hz. 
Subject headings: black hole physics — accretion, accretion disks — X-rays: binaries — X-rays: 
galaxies 



1. INTRODUCTION 

Non-nuclear ultraluminous X-ray sources (ULXs) in 
external g alaxies ha v e bee n known since the Ein- 
stein era fFabb iano '1989) with apparent isotropic 
luminosities above the Eddington limit of a IOMq 
black hole {L 2 x iQ^s erg/s) (for revie ws see 
IFahbiano Whitell2?M "Miller Cn]berdl200l. From 
the surveys of nea rby gal axie s und ertaken with the 
ROSAT HRI JColb ert & PtaF "200^) and the Chan- 
dra ACIS ijSwartz et al. .2004) , more than one hundred 
sources with luminosity greater than 10^^ erg/s have 
been located in nearby galaxies. 

The nature of ULXs is still unclear. The major dispute 
arises from our observational ignorance of the angular 
distribution of the source emission. If X-rays are emit- 
ted isotropically, then the high luminosity indicates that 
the ULXs might harbor intermediate mass black holes 
(IMBHs ) with the mass in the order of lOOMfn or even 
higher ( Colbert & Mushotzkv '1999'; 'Makishima ct al.' 
iOQQ; Ptak &_Griffiths 1999; Kaaret e t al. 2001). If the 
emission is beamed, then the ULXs could be stellar mass 
black holes ijKing et alJl200lt iKording et alJl2002) . 

Evidence that the luminosities are truly high for 
some ULXs has been found from their associations 
with optic al nebula which show indications of X-ray 
ionization ijPakull fc Mirionil [200^ . The photoioniz- 
ing flux required to power the nebula surrounding one 
ULX requires a lum inosity greater than 4 x 10'^^ erg/s 
l|Kaaret et al.ll2?)?)3) . Optical nebulae associated with 
some ULXs also show shoc k emission indicative of un- 
usually powerful explosions ijPakull fc Mirion'ill2003f) . 



X-ray spectral analysis has revealed thermal compo- 
nents in the s pectra of several ULXs with temperatures 
of 0.1-0.3 keV llColbert fc MushotzkvllT999t iKaaret et all 
120031 iMiller et al J 12003(1 . which are lower than that in 
stellar mass black holes in high states, suggesting accre- 
tion o nto a black hole with higher mass l)Makishima et all 
I2OOOI) . Also, a few ULXs exhibit breaks in their timing 
power spectra at relatively low frequencies which have 
been interpreted as evidence for high b lack hole masses 
l|CroPDer et all 120041 ISoria et alJl2004|) . However, such 
interpretations are ambiguous due to our lack of under- 
standing of timing power spectra and the fact that some 
stellar-mass black holes occasionally exhibit brea ks at 
similar frequencies l|McClintock fc Remillard]l2004|) . 

XMM is more powerful for spectral and timing analysis 
than previous X-ray telescopes owing to its large effective 
area. The goal of this paper is to use XMM archival data 
to study the spectral and timing behaviors of ULXs in 
nearby galaxies. In § El we describe the details of data 
reduction. We present the spectral analysis and details 
for some selected ULXs in § |31 and timing analysis and 
results in § 0] Finally, we interpret our results in § |S1 

2. DATA REDUCTION 

We derive our source sample from previously reported 
ULX dete ctions. Spec i fically , we used the Chandra 
catalog of 'Swartz et al.l l|2004|) . the ROSAT catalog of 
Colbert & Ptak ( 2002) , and the catalog of Liu & Mirabej 
(2005). We included only sources located in nearby 
galaxies with distances less than 11 Mpc. We selected 
sources with enough counts in XMM exposures for good 
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spectral and timing analysis. We eliminated sources 
crowded near other point sources that can not be resolved 
by XMM, e.g. the ULXs in M82. Starting from these 
three catalogs and applying the criteria above, which are 
solely intended to enable us to perform reliable spectral 
and timing analysis, we produced a final sample of 28 
sources, see Table ^ Hereafter we refer to sources by 
their number in Table^since not every ULXs has a com- 
mon name. Our sample is not meant to constitute an un- 
biased survey of nearby galaxies for ULXs. Rather, it is 
intended to provide a sample with which we can examine 
the spectral and timing properties of sources previously 
reported as ULXs in the literature. 

All data were processed using SAS 6.1.0 with the most 
recent calibration files as of March 2005. We extracted 
single and double events (PATTERN < 4) for pn data 
with the addition of triple events (PATTERN < 12) for 
MOS data. For spectral analysis, we required FLAG = 
and screened for good time intervals without background 
flares. For timing analysis, we required FLAG equal to 
#XMMEA_EM for MOS and #XMMEA_EP for pn, and 
selected a continuous good time interval without back- 
ground flare contamination or data gaps. We found that 
gaps in the XMM-Newton data stream often cause spu- 
rious timing signals in the frequency region 10~^ ~ 10^^ 
Hz. Therefore it is better to avoid gaps than attempt to 
fill them, e.g. with an average count rate. 

Source counts are generally extracted from a 32"- 
radius circle, following the recommendation in the 
data analysis threads on the XMM SAS web page at 
|http : / / xmm . vilspa .esa.es/ sas/ documentation/ threads/ 
For some very bright isolated sources, wc used a 40"- 
radius circle (this results in an increase in luminosity by 
a few percents but no difference in spectral parameters 
within errors). For sources in crowded regions or with 
a strong diffuse background, we used a source region 
with a 16"-radius circle. To attain the best statistics, 
we combine the data from the pn and the two MOS 
detectors unless some CCD detectors are turned off, the 
source rides on the gap between two CCD chips in one 
of the detectors, or the source is located at the edge of 
a CCD chip. The energy spectra are extracted from 
individual MOS and pn detectors and then grouped for 
fitting. The light curves are obtained from the merged 
events files. 

3. SPECTRAL ANALYSIS AND RESULTS 

The energy spectra are fitted with XSPEC 11.3.1. 
First, we fit all 28 sources with a single power-law 
spectrum and, separately, with a single multicolor disk 
(MCD) (Mitsuda ct al. 1984; Makishima ct al. 1981 
diskbb in XSPEC) spectrum. Both spectral models, 
and all others considered here, also include interstellar 
absorption. The absorption column density is treated 
as a free parameter with a lower bound set to the col- 
umn density along the line of sight within the Milky Way 
(jDickev fc Lockmanlll99d|) for all models in this paper. 
Then we tried adding a MCD component to each power- 
law model. We accept the MCD component if its addi- 
tion improves the fit by at least Aa significance as calcu- 
lated using the F-test for an additional component. The 
spectral parameters from all the fits, the unabsorbed 0.3- 
10 keV luminosities calculated assuming isotropic emis- 
sion, and x^/dof for each fit are listed in Tabled We 



include only those fits with x^/dof < 1.4. We find 13 
sources with a significant MCD component in addition 
to the the power-law component. In half of these sources, 
the power-law is dominant at high energies and the MCD 
disk temperature is in the range 0.1-0.4 keV. 

Source 1 is better fitted with a single MCD model than 
with a single power-law model. The F-test shows that 
the MCD model is preferred at a significant level of 0.05. 
For sources 2, 11, 13-16 and 19-21, the F-test shows no 
preference for a single MCD or a single power-law. For 
sources 9, 10, 18, 24, and 27, the x^/dof for the single 
MCD model is greater than 1.4, therefore that model is 
rejected, and no significant improvement was obtained 
with the MCD plus powerlaw model. For these sources, 
only the powerlaw fit is reported. 

The spectra of sources 12, 22 and 26 show significant 
residuals when fitted with a MCD plus power-law model 
and have x^/dof > 1.7. Therefore, we searched for an- 
other spectral model which would give a better fit. All 
three spectra show a bump around 0.5-1.0 keV, which 
is a region where hot diffuse gas produced significant 
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Fig. 1. — Energy spectra of source 22 fitted with an MCD model 
plus power-law model (top) or with a mekal plus power-law model 
(bottom). The mekal model is a significantly better to fit the 
spectrum because it adequately models the bump around 0.5-1.0 
keV. 
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Fig. 2. — Energy spectra of sources 17 and 7. The spectrum of source 17 (top left) consists of a 0.17 keV MCD and a power-law (top 
right). The spectrum of source 7 is curved at high energy (bottom left) which could be fitted either by a combination of a 1.51 keV MCD 
and a power-law (bottom right) or by a single Comptonization model (not shown). 



line emissi on. We found tha t the sum of a hot diffuse 
gas model (jMewe et al.lll995l mekal in XSPEC) plus a 
power-law produced a much better fit than the MCD plus 
power-law model for all three sources. A comparison of 
the two models for source 22 is shown in Fig. ^ 

For sources 3, 6, and 7, the MCD component is dom- 
inant at high energies. We fitted the spectra for these 
three sources with a single MCD m odel and with a sin- 
gle Comptonization model (Sunvacv fc Titarc huk 198(| 
compst in XSPEC). The single MCD model fits are not 
poor, with xi = 1-2 ^ 1-3, but are rejected compared 
to the MCD plus power-law model or the Comptoniza- 
tion model. The Comptonization model fits, with large 
optical depth r « 20, are as good or better than the 
power-law plus MCD model fits. These three sources 
shown significant curvature in the high energy parts of 
their spectra, which is quite distinct from the power-law 
followed by most of the other sources at high energies. 
Fig. 121 shows the spectrum of source 17, which is best 
fitted with an MCD plus power-law with the power-law 
dominant at high energies and which is typical of the 
sources with low disk temperatures, kT = 0.1 — 0.4 keV. 



In contrast, the spectrum of source 7 is the strongly 
curved at high energies, as shown fitted with an MCD 
plus power-law model in which the MCD component is 
dominant at high energies or a Comptonization model. 
The spectra of sources 3 and 6 have the same shape as 
source 7. 

In Fig. 131 we plot source luminosity L versus photon 
index F for sources with and without a significant ther- 
mal component. For sources without a thermal compo- 
nent there is no apparent correlation between luminosity 
and photon index. For sources with a thermal compo- 
nent, there appears to be an absence of sources with hard 
spectra at low luminosities, i.e. in the bottom left corner 
of the bottom panel of Fig. 13 

The relation between the photon index F of the power- 
law component and the temperature kT of the thermal 
component is shown in Fig. 01 for all sources where addi- 
tion of the thermal component significantly improves the 
fit. Considering all sources together, high temperatures 
appear to be correlated with soft photon indices. How- 
ever, this correlation appears to be a by-product of the 
juxtaposition of different spectral shapes. In particular. 
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Fig. 3. — Luminosity vs. photon index of ULXs. Top: ULXs 
with a single power-law spectrum; Bottom: ULXs with a prominent 
MCD component plus the power-law spectrum. The label indicates 
the source index number in Table IH 
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Fig. 4. — Photon index vs. temperature of ULXs with a promi- 
nent thermal component plus the power-law spectrum. The label 
indicates the source index number in Table [Tl 



sources 3, 6 and 7 which have the highest temperatures 
are well fitted by a single Comptonization model or a 
power-law plus MCD model with the power-law is domi- 
nant at low energies and the thermal component at high 
energies, in contrast to the other sources. The intermedi- 
ate temperature sources, 12, 22, and 26, are exactly those 
better fitted by a mekal model than a MCD model. Con- 
sidering only the sources for which the power-law plus 
MCD is the best fit and the power-law dominates at high 
energies, there is no correlation apparent between F and 
kT. 

Fig. |S1 shows the relationship between the temperature 
kT of the thermal component and the X-ray luminosity 
L, for all sources where addition of the thermal compo- 
nent significantly improves the fit. Again, the sources 



Fig. 5. — Luminosity vs. temperature of ULXs with a significant 
thermal component in their spectra. The label indicates the source 
index number in Table IH and the fractional Eddington luminosity 
in the brackets which is derived from the MCD model (for details, 
see the text). In the bottom panel, two ULXs which have been ob- 
served with XMM-Newton to undergo state transitions are added. 
The filled stars represent Holmberg II X-1 (Dewangan et al. 2004) 
and the white stars represent Holmberg IX X-f I Wang et al. 2004). 
The dashed lines connect individual sources observed in different 
states. 



segregate according to the best fitting spectral model. 
Sources 3, 6 and 7 have the highest temperatures and 
cover a wide luminosity range extending up to 10^° erg/ s. 
The sources with mekal components or a single MCD 
are at intermediate temperatures and relatively low lu- 
minosities. The sources for which the power-law plus 
MCD is the best fit tend to have low temperatures and 
relatively high luminosities. There is no correlation ap- 
parent between L and kT if we consider only the latter 
sources. 

We can, crudely, estimate the ratio of the observed 
luminosity to the Eddington luminosity using the fit- 
ted temperature and X-ray luminosity. By adopting the 
canonical values of t he color correction and te mperature 
profile correction in 'Makis hima et al] (|2p00), the disk 
blackbody temperature scales with black hole mass as 

fcT = 1.2r^.^(^j keV, (1) 

where /3 is determined by the black hole spin, (3 = 1 for 
a Schwarzschild black hole and /3 = 1 /6 for a maximally 
rotating Kerr black hole, and 7] is the bolometric lumi- 
nosity normalized to the Eddington luminosity as 77 = 
iboi/iEdd- We take Ledd = (1-3 x lO^Serg/s) (M/Mq) 
and approximate Lboi = -^x- Solving for 77, we find the 
fractional Eddington luminosity 
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We derive a fractional Eddington luminosity (shown 
in Fig. in the brackets) from Eq. 13 with (3 = 1. There 
are three sources which appear to be above the Edding- 
ton hmit, with L/L^dd of 1.2, 1.7 and 4.1 respectively 
for sources 3, 6, and 7. If the black hole s are r apidly 
rotating, as suggested bv iMakishima et al.l (|2000), then 
the sources fall below Eddington luminosity with source 
7 still near the Eddington luminosity. 

4. TIMING ANALYSIS AND RESULTS 

After examining the power spectral density (PSD) of 
every source in our sample, we find three sources (No. 4, 
23, & 25) with detectable noise power in the frequency 
range 0.001-0.1 Hz. The lower frequency bound is chosen 
based on the shortest good time interval of any observa- 
tion used in this analysis in order to allow a uniform 
frequency range for all sources. The PSD in the other 
sources is consistent with Poisson fluctuations, and we 
can only place an upper limit on the noise power (see 
Fig. EJ. There is no obvious correlation between the rms 
variability and the luminosity among those three sources 
with detectable noise power. The rms upper limits ap- 
pear correlated with source luminosity because the sen- 
sitivity of the upper limit depends on the total number 
of counts recorded. 

We are able to detect timing noise in several frequency 
bands for three sources. No. 7, 22, & 25. Light curves 
and PSDs for these three sources are presented in F ig. [7| 
All PSDs are normalized to rms ijvan de KHsllTosaj) and 
binned in logarithmic scale with points above the Poisson 
level plotted with error bars and points below the Poisson 
level plotted as an upper limit. We have checked the PSD 
from background regions for these observations to ensure 
that the timing noise is not produced by variations in the 
background. The noise power for sources 7 and 22 lies 
below 0.001 Hz. They do not have detectable noise power 
in 0.001-0.1 Hz band used above. 

We note t hat the timing features reported by 

ICropper et al.l l)2004|) . including PSDs with a break 



power-law in source 17 (NGC 4559 X7) and a single 
power-law in source 18 (NGC 4559 XIO), are not re- 
peated in our results. We examined the same XMM 
observation and found multiple background flares dis- 
tributed across the whole exposure. If we only extract 
light curves at clean intervals, even including several 
small flares and with a duration up to 80% of the whole 
exposure, no timing noise is present above the Poisson 
level. If we use the whole exposure including the back- 
ground flares, we ca n reproduce the PSD reported by 
iCropper et al.l l)2004|) . We con clude that the f eature s in 
the power spectra reported bv lCropner et al.1 l)2004l) are 
caused by background flares. 

No. 7 (NGC 3031) We calculated the PSD from a 
40.6 ks good time interval without timing gaps or back- 
ground flares for both the M0S2 and the pn CCDs in 
the first half of the observation. The second half of the 
observation contains strong background flares and was 
discarded. The PSD is calculated in the energy range 2- 
10 keV with a 4096-point FFT. This energy band was 
selected because the timing noise is strongest in this 
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Fig. 6. — The relationship between the fractional rms integrated 
in the frequency region 0.001-0.1 Hz and the X-ray luminosity 
(top), power-law photon index (middle) and the thermal tempera- 
ture (bottom). The luminosity is taken from the best-fitted spec- 
trum. Three sources have detectable timing noise; these sources 
are best fitted with a MCD plus power-law spectrum. The sources 
for which timing noise above the Poisson level is not detected are 
shown as 2-cr upper limits. In the middle and bottom panels, only 
sources which are well fitted with a MCD plus power-law model 
are plotted. Sources without detected noise are plotted only if the 
upper limit on the rms is less than 0.05. The diamonds indicate 
sources which were well fitted with a Comptonization model or a 
high temperature MCD plus power-law model. 



band. For sources 22 and 25, we also present the PSD 
for the energy which maximizes the timing noise. The 
time resolution of the light curve is set to the interval 
duration divided by the number of FFT points (as for all 
these three sources). Power at frequencies below ~ 0.1 
mHz increases at lower frequencies with power-law slope 
a = -1.8 ± 0.3. The PSD below 2 keV only has the 
lowest frequency point above the Poisson level. 

No. 22 (NGC 5128) We select a 11.4 ks good time 
interval without gaps or background contamination for 
three MOS and pn CCDs. The PSD is calculated in 
the energy range 0.7-10 keV with a 1024-point FFT. We 
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Fig. 7. — Light curves and power spectra of three ULXs. 
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also examined a power spectrum for events with energies 
below 0.7 keV, but there is no detectable power. The 
power spectrum below 0.2 mHz has a power-law slope 
a = — 1.56±0.14. There is also some power at frequencies 
around 2 mHz, but our sensitivity is inadequate to make 
unambiguous statements about this. 

No. 25 (NGC 5408) The observation contains 
strong background flares in the second half and a 2.96 
ks good time interval in th e first half, w hich is used for 
calculation of the PSD. .Soria et alJ l|2004D report a power 
spectrum for this observation with a power-law at high 
frequency breaking to a flat continuum (a = 0) below a 
cutoff frequency of 2.5 mHz. Such a PSD is often found 
in Galactic black hole binaries and AGNs and the break 
frequency scales with mass, although the break frequency 
in individual sources varies by factors of at least 10. We 
re-examine the PSD in the 0.2-10 keV band with a 1024- 
point FFT. At frequencies below 0.001 Hz the power 
spectrum has a power-law slope a — —0.8 ± 0.2. Above 
this frequency, there appears to be power distributed ir- 



regularly across a range assoc i ated w ith the flaring time 
scale reported bv 'Soria et alJ l)2004f) . Our power spec- 
trum does not suggest a break to a flat continuum with 
a = at low frequen cy. We are able t o reproduce the 
spectrum reported bv lSoria et al.l l)2004|) if we use much 
broader frequency bins. 

5. DISCUSSION 

We have found a diversity of spectra in our sample 
of ULXs. Considering only those sources for which the 
statistics justify the use of models more complex than a 
simple power-law, half the sources are best fitted with 
a MCD plus power-law model with the power-law dom- 
inant at high energies. This model is commonly used 
in fitting the spectra of accreting black holes. These 
sources have the lowest thermal component tempera- 
tures, 0.1-0.4 keV, and extend to the highest luminosi- 
ties, above 10^° erg/s. Three sources have spectra which 
are strongly curved at high energies. These sources are 
best fitted with a MCD plus power-law model with the 
MCD component dominant at high energies and a steep 
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power-law at low energies or a Comptonization model. 
These sources have the highest temperatures in our sam- 
ple, 1.0-1.4 keV. Three sources show evidence of diffuse 
hot plasma and are best fitted with a mekal model plus 
a power-law. The temperatures of these three sources 
are in the range 0.6-0.8 keV and the luminosities are the 
lowest in our sample, near 10'^^ erg/s. Finally, there is 
a single source best fitted with an MCD model alone. 
This source is likely a stellar-mass black hole in a high 
state. Within each class, the spectral parameters, par- 
ticularly the emission temperature of the thermal com- 
ponent, appear to be rather similar. Between classes, 
the temperature of the thermal component varies signifi- 
cantly. This diversity of spectral properties and the fact 
that the sources appear to lie in three distinct classes, 
see Fig.O suggests that the ULXs are, in fact, a diverse 
population. 

The three apparently distinct classes of ULXs may 
merely represent three states of ULXs and not a true 
diversity of physical objects. To examine if this is the 
case, we searched the literature for state transitions of 
ULXs in which the temperature of the thermal com- 
ponent has been measured . We fo und four such cases: 
Holm berg II X-1 ijP ewang an et al]l2 004). Holmberg IX 
X-1 La Parola et al. (2001); Wang et al. m2004|) . and two 
ULXs in IC 342 ( Kubota et al. 2001). For Holmberg II 
and Holmberg IX X-1, multiple observations with XMM- 
Newton are available. We plot the spectral states of these 
sources, connected by dashed lines, in a bottom panel of 
Fig. |S1 As can be seen from inspection of the figure, the 
ULXs which undergo state transitions remain within the 
same group in the luminosity vs temperature diagram. 
The state transitions do not move the sources between 
the groups. Both sources remain in the low kT through- 
out the state transitions. For the state transitions of the 
two ULXs in IC 342 reported by Kubota ct al. ( 2001), 
both sources remain in the high kT group throughout the 
transitions. These sources are not plotted in Fig. [S] be- 
cause we prefer to plot only spectra from XMM-Newton 
observations. 

We note that lLa Parola et al.l lj200H) reported the state 
transition of M81 X-9 = Ho IX X-1 (source 8) with a 
high state observed by ASCA fitting by a single MCD 
model with a high temperature kT — 1.24 keV and a 
low state fitted by a power-law plus MCD model with 
kT ~ 0.3 keV. If correct, this would represent a tran- 
sition between__the low A;T and high kT groups. How- 
ever, IWang et all ll200l reported two observations us- 
ing XMM-Newton of Ho IX X-1 (one is the same as 
used in this paper). The luminosity during the XMM- 
Newton observations is a factor of 2 higher than in the 
ASCA observati on and clear l y repr esents the high state 
of the source. iWang et alJ ll2004D concluded a single 
MCD model, as used bv ILa Parola et alJ l|2001|) to fit 
the ASCA data, is strongly ruled out by the high qual- 
ity XMM-Newton data. Instead, the spectra are better 
described by a cool disk plus power-law model. Their 
results are consistent with our own which rule out appli- 
cation of a single MCD model to the high state of Ho IX 
X-1. Therefore, we conclude that Ho IX X-1 remains in 
the low kT in its high and low states. 

Those ULXs which undergo state transitions remain 
in the same class throughout the transition, suggesting 
that the definition of the different classes is robust. We 



note that, on average, the low temperature ULXs are 
more luminous than the high temperature ones. There- 
fore, if the different classes represent different states of 
a single underlying population, then the lower luminos- 
ity states would be associated with higher temperatures, 
which is the reverse of the normal situation for black hole 
X-ray binaries. However, additional XMM-Newton ob- 
servations of the sources in our sample would be of great 
interest to test if these source classes truly represent dis- 
tinct classes of objects. 

5.1. MCD plus power-law sources 

Here, we consider the properties of the class of sources 
which are best fitted with MCD plus power-law model 
with the power-law dominant at high energies. In this 
class, there is an absence of sources with hard spectra 
at luminosities near 10^^ erg/s. This is not the case 
for sources which lack a detectable thermal component. 
When a black h ole is in its low/hard state, the disk tem- 
perature is low (fM cClintock fc R emillard 2004) . If these 
ULXs are intermediate mass black holes (IMBHs), then 
their temperature in low/hard states would shift out of 
the detection threshold of XMM. This would naturally 
produce an absence of sources with hard spectra at (rel- 
atively) low luminosities. This scenario would imply that 
the sources with single power-law spectra with hard spec- 
tra (r < 2) and luminosities near lO'^^ erg/s may be 
IMBHs in the low/hard state. 

For accreting black holes in the high state, higher 
black hole masses should lead to cooler disk temperatures 
l^akishima ct al. 2000). The low temperatures of some 
ULXs have b een suggested to be evidence t hat they con- 
tain IMB Hs llColbert fc Mushotzkvl [19991 iKaaret et all 
[2003; Mill er et al.ll2003(l . In Fig. 151 the sources which are 
best fitted with a MCD plus power-law model with the 
power-law dominant at high e nergies concentrate at low 
temperatures, kT < 0.4 keV. IMiller et all l|2004n have 
suggested that the clustering of ULXs in a narrow range 
in temperature is evidence that those ULXs are IMBHs 
with masses lying in a narrow range. However, current X- 
ray instrument lack the sensitivity to detect thermal com- 
ponents with temperatures much below 0.1 keV. Hence, 
the lower boundary of the range is observational. The 
thermal emission from higher mass black holes would 
simply not be detected. 

In stellar mass black holes, different spectral states are 
associated with different levels of variability. The to- 
tal noise power integrated in the 0.1-10 Hz band has 
an rms < 0.06 for high states and an rms of 0.1 to 0.3 
for low states (jMcClintock fc Remillardll200^ . Because 
these relatively high frequencies are not accessible for 
ULXs given the sensitivity of XMM-Newton, we use the 
0.001-0.1 Hz band to calculate the rms noise power. The 
relations between the rms noise level and the luminosity, 
thermal component temperature, and power-law photon 
index are shown in Fig. El For those three sources with 
detected rms, the noise does not appear correlated with 
luminosity. Also, the noise level does not appear corre- 
lated with photon index. Indeed, the one source which 
has a high noise level, as associated with the hard state of 
stellar-mass X-ray binaries, has a relatively soft photon 
index. 

The X-ray variability of Galacti c X-ray binaries at 
low frequencies has been examined bv lGilfanov fc ArefievI 
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()2005tl . In the frequency region from 10^"' to 0.1 Hz, the 
PSD of X-ray binaries is found to follow a power-law with 
an exponent a in the range —1.2 to —1.5. The PSDs for 
source 7, 22 and 25 at low frequencies follow power-law 
forms with exponents a of —1.8 ±0.3, —1.56 ±0.14, and 
0.8 ± 0.2, respectively. These are similar to the PSDs 
of Galactic X-ray binaries at low frequencies. This may 
suggest that the PSDs reflect viscous variabilities from 
the accretion disks of the ULXs. 

5.2. Mekal sources 

The sources 12, 22 and 26 are not adequately fit- 
ted with a MCD plus power-law model. Their spectra 
present a prominent bump around 0.5-1.0 keV and their 
spectra are adequately fitted with the sum of a mekal 
thermal plasma model plus a power-law. The presence 
of a mekal component suggests that part of the X-ray 
emission from these sources arises from a hot, diffuse 
thermal plasem. The fractions of the total luminosity 
in the mekal component are 9.3%, 32.2% and 31.7% for 
source 12, 22 and 26, respectively. 

We can place an upper bound on the physical size of 
these sources from that fact that all three sources ap- 
pear unreso l ved w ith Chandra: Obs ID 402 for source 12 , 
iKraft et alJ l)2001|) for source 22, and lSwartz et all 1)20041) 
for source 26. Therefore their angular size must be less 
than 0'.'5. With this bound on the physical size and the 
mekal fit parameters temperature kT and normalization 
K , we can then derive a lower limit on the electron den- 
sity n > ^J'iTTD'^K/{10^'^V) (cm~^), an upper hmit on 
the total energy of i? < 2 x ^kTnV (erg), and an upper 

limit on the cooling time oi t < 1.05 x lO^^/kT/n (yr), 
where D is the distance to the host galaxy in cm, V is 
upper limit of the plasma volume V < It"^ in cm^, kT is 
the plasma temperature in keV, and the value and units 
of K are taken from XSPEC. 

No. 12 (NG C 4395) T he source is also known as 
NGC 4395 X-2 l|Lira et all 2000) located in the dwarf 
galaxy NGC 4395. In the XMM observation on 2003- 
Nov-30, the source is located near the gap in pn CCDs, 
leading to a small difference in flux between pn and MOS 
CCDs. Therefore we only take the two MOS CCDs for 
spectral analysis for this source. The spectral fitting 
gives kT = 0.78 ± 0.04 keV and = 4.5 ± 0.5 x 10-^ 
from which we derive n > 9.3 cm^'^, E < 2.8 x 10^"*^ erg 
and t < 1.0 X 10^ yr. 

No. 22 (NGC 5128) The source is located near 
the edge of southwest radio lobe in NGC 5128 (Cen- 
taurus A). Our result shows the source has a strong 
mekal component. The derived plasma properties from 
the temperature kT = 0.68 ± 0.03 keV and normal- 
ization ivT = 4.1 ± 0.04 X 10"^ are n > 7.6 cm~^, 
£; < 5.1 X 10^1 erg and t < l.l x 10^ yr. The X-ray 
source, according to the coordinates RA=13h25m 07.45s, 
Dec= -43d04m09.53s (J2000) from Chandra (Kra ft et all 
1200 1|) . corresponds to an optical co unterpart that is in- 
cluded in the USNO-Bl.O catalog l|Monet et alJ 1200^ 
with an apparent magnitude mB2 = 14.23 mag and 
'7ifl,2 ~ 14: Al mag, corresponding to an absolute mag- 
nitude Mb2 = -14.2 mag and Mji2 = -14.0. The B-R 
color of the source is i nconsistent with a foreground M 
dwarf as suggested by l|Kraft et al.ll200jl) . 

No. 26 (NGC 5457) The source is located 6' west of 



the nucleus of NGC 5457 (MlOl), with a relatively low 
luminosity {L < 10^^ erg/s for the best-fitting model) 
in this XMM observation. The mekal component is 
prominent in this observation with spectral parameters 
kT = 0.67±0.04 keV and K = 1.7±0.2 x 10"^. We derive 
the plasma parameters n > 4.0 cm~^, E < 7.8 x 10^^ erg 
and t < 2.1 x 10^ yr. This source also has an optical 
counterpart in the DS S image. The optical counterpart 
is in the GSC catalog i Morris on et al. 2001) and has an 
apparent magnitude my = 15.5±0.3 mag, corresponding 
to an absolute magnitude My = —13.7 mag. This source 
is also observed by the Optical Monitor on-aboard XMM 
with an absolute magnitude Myvwi = —10.8 mag. 

One plausible interpretation for the hot diffuse plasma 
components is that they are young supernova remnants. 
The derived density, total energy, and cooling time are 
consistent with this interpretation. The X-ray luminosi- 
ties and temperatures of these three s ources are similar 
to rec ent supernovae (e.g . SN 1979C, llmmler & LewiiJ 
120031 llmmler et aLll200,5fi . The high absolute magni- 
tudes of the optical counterparts indicates they are not 
normal stars. Comparison with the late-time absolute 
magnitude of SN 1979C (Mb = -12.1 mag, iVfuvwi — 
— 13.6 mag), suggests that they may be supernova rem- 
nants which are several tens years old. We do not find 
supernovae counterparts on the lAU supernova list^, but 
the supernova may easily have been missed. Source 22 is 
found to have noise power above the Poisson level, which 
could not be produced by a spatially extended source. 
However, more than 2/3 of the luminosity arises from the 
power-law emission. The low frequency noise for source 
22 (the first two points in Fig.[7I) is strongest in the 1-10 
keV band, where the energy spectrum is dominated by 
the power-law component also. Therefore, the variability 
appears related to the power-law rather than the mekal 
component and may arise from a central compact object. 

5.3. High temperature sources 

Three sources (3, 6, and 7) have spectra which are 
strongly curved at high energies in strong contrast to 
most of the other sources. From the MCD fits, we found 
that the sources have high temperatures and are either 
super-Eddington for a stationary black hole or near Ed- 
dington for a rapidly rotating black hole. These unusual 
aspects establish these sources as a distinct class. All 
three sources have low timing noise with an rms in the 
0. 001-0.1 Hz band belo w 0.02. 

iFoschini et'all l)2004|) suggested that source 3 (M33 X- 
8) is a stellar-mass black hole in the high state because 
its behavior is similar to that of Galactic black hole X- 
ray binaries. Radio emission with a spectral index simi- 
lar to microquasars h as been found in VLA observations 
ijGordon et alJ |l ^^^. po ssiblv linking this source with 
microquasars. iFoschinie t al. (2004) also suggest that the 
emission may be beamed. Source 7 (in M81, F88 X-6) 
has also been identified with a we ak radio counterpart 
(|Fa,bbia,nnll198aiSwaTtz et a,lj[200l . 

Other interpretations for the high temperatures in 
some ULXs i nclude intermediate mas s black holes in the 
Kerr metric (|Makishima et all 120001) or black holes of 
several tens of solar masses with slim accretion disks 
(advection dominated optically thick disks) emitting 

^ |http : //cf a-www ■ harvard ■ edu/iau/list s/Supernovae ■ html | 
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at super-Eddington luminosities ijEbisawa et al.l l2003j) . 
These interpretations are based on results showing that 
single multicolor disk blackbody spectrum provi de the 
best fits for high-temperature ULXs (Makishim a et all 
l200n|) . Our results based on the high quality spectra 
available from XMM-Newton show that a single MCD 
provides a significantly worse fit than a Comptoniza- 
tion model. A MCD model with an additional power- 
law component can provide an equally good fit, but the 
power-law is dominant at low energies, which is differ- 
ent from the behavior usually seen from Galactic black 
hole binaries where the power-law component is usu- 
ally weak when then source is in the thermal dominant 
state which shows prominent M CD emission near 1 keV 
ijMcClintock fc RemillardI l2004j) . Interpretation of the 
low-energy part of power-law as due to Comptonization 
would be problematic because the seed photons upscat- 
tered by Comptonization are usually assum ed to arise 
from the disk emission (Ro berts et al.ll2'005fl . Power-law 
components extending below the disk emission have been 
found in the spec tra of the b lack hole candidates LMC X- 
1 and LMC X-3 ijHaardt et al. 2001), but are likely bet- 
ter interpreted in terms of physically motivated Comp- 
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tonized MCD models l)Yao. Wang fc Zhangl2005j) . In the 
Comptonization model, a high electron column density 
~ 3 X 10^^ cm~^, is req uired to produ ce the observed op- 
tical depths l|Makishim a et all I200CI() . Since heavy ions 
would not be ionized at the inferred temperatures, strong 
Fe K absorption should be observed unless the material 
has very low metallicity. No Fe K absorption is appar- 
ent in the data, but the constraints are weak due to the 
limited statistics at high energies. Adding a smeared 
edge model to the Comptonization spectra with a cen- 
tral energy of 7 keV, optical depths up to r ~ 6 are 
allowed, depending on the edge width. Lack of Fe K ab- 
sorption would argue against accepting the Comptoniza- 
tion model. 
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TABLE 1 

XMM Observations of ULXs in nearby galaxies 



No. Host Galaxy 


RA 


Dec 


Dist. 


Obs Date Exposure'' Ref.'^ 




(J2000) 


(J2000) 


(Mpc) 


(ks) 



I 


Circinus 


14 


12 


53. 


.52 


—65 


22 


54, 


7 


4.0 


08/06/01 


61.1/91.1 


I 


2 


IC 342 


03 


45 


55. 


7 


+68 


04 


58 




3.9 


02/1 1 /01 


4.9 /5.1 


2 


3 


NGC 


598 


01 


33 


50. 


.89 


+30 


39 


37, 


2 


0.9 


08/04/00 


6.5 /II. 9 


I 


4 


NGC 


1313 


03 


18 


19. 


,9 


—66 


29 


07 




3.7 


10/17/00 


17 4/28 n 


2 


5 


NGC 


1313 


03 


18 


22. 


2 


—66 


36 


01 




3.7 


10/17 /OO 


1 7 4/28 n 


2 


6 


NGC 


2403 


07 


36 


25. 


.55 


+65 


35 


40, 





3.1 


09/12/04 


50 6/69 1 


1 


7 


NGC 


3031 


09 


55 


32. 


.98 


+69 


00 


33, 


4 


3.6 


04/22/01 


87 8/125 4 


I 


8 


NGC 


3031 


09 


57 


54. 


,4 


+69 


03 


43 




3.6 


04/16/02 


7.7/7.7 


2 


9 


NGC 


3628 


11 


20 


15. 


.76 


-H 3 


35 


13, 


7 


10.0 


1 1 /27 /on 


37.7/41.2 


I 


10 


NGC 


4258 


12 


18 


43. 


.86 


+47 


17 


31, 


,5 


7.7 


1 2/08/00 


12.2/14.4 


1 


11 


NGC 


4258 


12 


18 


57. 


.84 


-1-47 


16 


07, 


I 


7.7 


1 2 /OS /on 


12.2/14.4 


I 


12 


NGC 


4395 


12 


26 


01. 


,8 


+33 


31 


34 




3.6 


1 1 /30/03 


86 9 /96 1 


2 


13 


NGC 


4449 


12 


28 


17. 


.83 


-1-44 


06 


33, 


9 


3.7 


06/02/02 


10 fi/1 1 8 


X 


14 


NGC 


4485 


12 


30 


30. 


.56 


+41 


41 


42, 


,3 


7.8 


05/27/02 


11.5/12.0 


1 


15 


NGC 


4490 


12 


30 


36. 


.32 


+41 


38 


37, 


,8 


7.8 


05/27/02 


11.5/12.0 


1 


16 


NGC 


4490 


12 


30 


43. 


,25 


+41 


38 


18, 


,4 


7.8 


05/27/02 


11.5/12.0 


1 


17 


NGC 


4559 


12 


35 


51. 


,71 


+27 


56 


04, 


,1 


10.3 


05/27/03 


33.7/36.4 


1 


18 


NGC 


4559 


12 


35 


58. 


,57 


+27 


57 


41, 


,8 


10.3 


05/27/03 


33.7/36.4 


1 


19 


NGC 


4631 


12 


41 


55. 


,56 


+32 


32 


16, 


,9 


7.6 


06/28/02 


44.9/53.2 


1 


20 


NGC 


4945 


13 


05 


21. 


,94 


-49 


28 


26, 


,6 


3.7 


01/21/01 


17.9/17.9 


1 


21 


NGC 


4945 


13 


05 


32. 


,89 


-49 


27 


34, 


1 


3.7 


01/21/01 


17.9/17.9 


2 


22 


NGC 


5128 


13 


25 


07. 


,4 


-43 


04 


06 




4.9 


02/02/01 


16.8/16.8 


2 


23 


NGC 


5204 


13 


29 


38. 


,61 


+58 


25 


05, 


,6 


4.5 


01/06/03 


15.4/15.4 


1 


24 


NGC 


5236 


13 


37 


20. 


,1 


-29 


53 


46 




4.7 


01/27/03 


18.9/24.1 


2 


25 


NGC 


5408 


14 


03 


19. 


,61 


-41 


22 


59, 


,6 


4.8 


01/28/03 


2.7/4.0 


3 


26 


NGC 


5457 


14 


02 


29. 


,89 


+54 


21 


19, 


,0 


7.0 


06/04/02 


21.6/36.5 


1 


27 


NGC 


5457 


14 


04 


14. 


,28 


+54 


26 


03, 


,6 


7.0 


06/04/02 


21.6/36.5 


1 


28 


UGC 


4305 


08 


19 


28. 


,99 


+70 


42 


19, 


,4 


3.4 


04/10/02 


4.7/4.7 


1 



''The clean/total exposure of pn CCDs, or mos CCDs in case pn data are not used. 
''The catalog from which the coordinates of the ULXs and the distance of the ho st gal axy 
are taken; fl') 4Swartz eTaTI <200l) . (2')- 4Colbert fc PtakI t2003) . f3') 4Liu fc MirabeB ^2005^ . 
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TABLE 2 

Spectral parameters and variabilities of selected ULXs. 



No. 


Model 


/ a 


JlH 


/ b 
1 /T 




r d 


2 /J f 
X /dof 


rms"^ 






l^iu cm f 


l^iu cm ) 




^Kev j 


(iO erg/s; 




t^U.UUl— U.i rLz) 


1 


powerlaw 


5.33 








42+0-01 
^^-^^-0.14 


340.4/265 


< 0.131 




diskbb 




5 fi+0.4 


^•"'-0.05 


74+0.03 
'^•'^-0.03 


0764+0-0002 

U.U/Dt_Q QQg7 


267.4/265 




2 


powerlaw 


3.03 


t- O+0.6 
3-«-0.6 




44+0-04 
'-'■^^-0.05 


142.9/153 


< 0.069 




diskbb 




o 1+0.4 
■^•-^-0.1 




^■^1-0.15 


n 4S+0 01 


169.2/153 




3 


diskbb 


0.56 


"^6+0 02 
'^•'^°+0.00 

1 7+0.2 
^•'-0.2 
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156+0-001 

U.10D_Q QQ3 
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< 0.003 




powerlaw+diskbb 




9 34+0.16 
^•■^^-0.13 


1 14+O.O6 
1-1^-0.06 


226+0-006 
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1 6+°l 




1 -,0+0.05 
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4 


powerlaw 


0.39 




1 05+0.04 




59+0-02 
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0.028 + 0.013 




powerlaw+diskbb 




'^•'-'-0.3 


1 74+0.06 
^•'*-0.06 


20+0-02 
'-'■^'-'-0.02 


n 01 +0.03 


598.9/625 




5 


powerlaw 


0.39 
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^•'-0.2 
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^•^'^-O.OB 




29+0-01 
'J-^^-o.oi 


266.9/252 


< 0.077 




powerlaw+diskbb 




2 9+°-^ 


2 

^•-■-^-0.15 


r, 90 + 0.05 

'-'■^■-'-0.06 


sn+0-02 

'^•'''^-0.05 


244.8/250 




6 


diskbb 


0.41 






1 04+0-02 
l-'J^-0.02 


170+0.001 
'^•l'^-0.003 
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''Absorption column density of the Milky Way IDickev Hi Lockmanlll990f) 

means the photon index for models including a power-law component; r means the optical depth for the compst model. 
'^kT stands for the inner disk temperature for the diskbb model, the plasma temperature for the mekal model, and the electron 
temperature for the compst model. 
'lUnabsorbed X-ray luminosity in the energy range 0.3-10 keV. 

•^rms variabilities in the frequency range 0.001-0.1 Hz. The rms values are not associated with the spectral models, while in 
Fig. |H]the luminosity is selected from the best-fitted model. 



